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We study cosmological solutions in the context of 4-diinensional low energy Heterotic M-theory 
with moving bulk branes. In particular we present non-trivial, analytic axion solutions generated 
by new symmetries of the full potential-free action which are similar to 'triple axion' solutions 
found in Pre-Big-Bang (PBB) cosmologies. We also consider the presence of a non-perturbative 
superpotential, for which we find cosmological solutions with and without a background perfect 
fluid. In the absence of a fluid the dilaton and the T-moduIus go to the potential-free solutions at 
late time, while the moving brane tries to avoid colliding with the boundary and stabilize within the 
bulk. When the fluid is included, we find that the real parts of the fields track its behaviour and that 
the moving brane gets stabilized at the middle point between the boundaries. In this latter case we 
^^ ' can make analytic approximations for the evolution of the fields whether or not a^xions are included 

f^ ' and we consider the possibility of this set up being a realization of the quintessential scenario. 

o' 

^ , PACS numbers: ll.25.Mj, 11.25.Yb, 98.80.Cq 



(N 



I. INTRODUCTION 



In recent years, there has been growing interest in the role played by extra dimensions in Early Universe cosmology 
In particular, a lot of work has been devoted to studying the so called brane-world cosmology, where the Universe is a 
4-dimensional extended object within a higher-dimensional space [IIQjBIJI- One particular, interesting, possibility is 
C , to consider moving branes 0, B S S S HO^ HB IH? \2M ■ I^ this paper we will look at moving bulk brane cosmologies 
^~^ • within a well defined M-theory context. 

jSt I To be more explicit, our starting point is the 5-dimensional effective action that has been derived by compactifying 

^j ' 11-dimcnsional heterotic M-theory [l4lll5l | on a Calabi-Yau 3-fold [ll lla IitI llM Il9| . The resulting theory is an explicit 

("^ [ realisation of a brane world scenario and consists of a 5-dimensional N=l su per grav ity theory on the orbifold S^/Zo 

coupled to two N=l 4-dimensional theories at the orbifold fixed points [H [lE ISS EllS • It has been shown [ilia. [la. 

[1^123,123 that one is free to add M-theory 5 branes into the 11-dimensional theory. Upon compactification, these 

—^ , branes appear as 3-branes in the 5-dimensional bulk and are free to move along the orbifold direction. To simplify 

i^ ' the problem, we will use the 4-dimensional effective action obtained [23. l23j by reduction on the 5-dimensional BPS 

r^ , domain wall vacuum solution. 

r* ' The minimal version of this set up contains three chiral superfields, namely the dilaton S, the universal T modulus 

. !^ I and the modulus Z, which specify the position of the moving brane in the orbifold direction. These superfields 

^^ . are constructed from six real scalar fields, three of which have a geometrical interpretation in terms of the higher 

i—i ' dimensional theory. 

. . . 1 This scenario has been studied, in the absence of a potential, in Ref. pdl |. where, it was shown that the axionic 

components of the chiral superfields could be set to zero. Then two of the remaining component fields, ip (measuring 

the Calabi-Yau volume) and the real part of the T modulus, /3 (measuring the orbifold size), behave asymptotically 

as rolling radii (RR) fields [2J|, while the bulk brane, z, moves a finite distance to generate a transition between 

the early and late time rolling radii solutions. These solutions are familiar from Pre-Big-Bang (PBB) cosmology and 

provide a framework for negative time branch cosmology within Heterotic M-Theory. 

In this paper, we will go further and take into account both new symmetries of the potential-free action, in order to 
generate new axion-dependent solutions, and also non-perturbative effects |2a l2a 1271 l28l l29l l30l l3ll| . In particular, 
we will introduce a superpotential due to M2 instantons stretching between the boundaries and the moving brane, 
and study the effects of the scalar potential coming from this superpotential in the evolution of the moduli. 

The outline of the paper is as follows. In the next section, we will explain our scenario in more detail, motivate the 
introduction of this particular superpotential and present the action we will work with and its symmetries. 
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In section IIIII wc will show how, in the absence of a potential, these symmetries generate new solutions, with 
non-trivial axion dynamics, from the ones found in Ref. |ll| . These solutions are valid on either time branch but 
might be particularly important in generating scale invariant isocurvaturc perturbations in a similar fashion to the 
PBB scenario. 

In section Hvl we see that we can still consistently truncate off the axions and investigate the behaviour of ip, j3 and 
z in the presence of the membrane potential on the negative time branch. We find that, at late times, the first two 
fields exhibit a subset of the late time rolling radii behaviour of the potential free case, but with the potential acting 
to prevent a bulk-boundary branc collision and possibly providing a way of transferring isocurvaturc perturbations 
into adiabatic ones (analogously to what happens, for example, in Ref. |32|'). 

We will, in section |V1 introduce a background fluid into the system, as required for a realistic model of the Universe 
evolving in the positive time branch. Firstly we will study its effect on the real parts of the moduli by setting the 
axions to their vacuum values. We will sec how the flelds tend to scale with the background fluid, allowing us to find 
a set of analytic solutions in different limiting cases. We also consider the possibility of one of these moduli being 
the quintessence field that might be accelerating the Universe today. Finally, we will include the axions' evolution 
towards their minima in the analysis and we will show how this will not affect the late scaling behaviour of the real 
parts of the moduli. 

We present our conclusions and future work in section IVll 

II. SET UP 

Our 5-dimensional branc- world scenario comes from the compactification on a Calabi-Yau three fold [l|, [Iflj llTl Ha . 
[l^|22| of the low-energy limit of 11-dimensional Horava-Witten theory [l4Lll5l| that is, 11-dimensional supergravity 
on the orbifold S^ IZ2 x Mio coupled to two 10-dimensional E% Super- Yang-Mills theories, one residing on each of the 
two 10-dimcnsional orbifold fixed planes. This compactification leads [iL Il9ll20ll2lll22 | in the bulk, to a 5 dimensional 
N=l gauged supergravity on the orbifold S^ /Z2 x M4 coupled to two N=l gauged theories on the now 4-dimensional 
orbifold fixed planes on which two 3-branes carry the so called observable and hidden sectors. 

In addition to this minimal set-up, we are free to include M-theory 5-branes in the vacuum of the 11-dimensional 
theory [l9l l23j . Upon compactification two of their dimensions wrap holomorphic two cycles within the Calabi-Yau 
three-fold, leaving effective 3-brancs in the 5 dimensional bulk. These additional 3-branes are transverse to the orbifold 
direction and have the interesting property of being able to move along that direction. Each of them will carry an 
additional N=l supcrsymmetric theory. For the purpose of our paper, wc will restrict ourselves to the inclusion of 
only one additional 3-brane. The charges on the visible and hidden orbifold planes and the three brane, respectively 
€[1 (i=l,2,3), are quantised in units of eg a-nd will satisfy the cohomology condition [33J 



E- 
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which follows from anomaly cancellation in the 11-dimensional theory. For the compactification to preserve 4- 
dimensional supersymmetry, g = (73 > [221 . 

Practically, we will be working in the context of a 4-dimensional effective theory which arises from further reduction 
of the 5-dimensional theory on the domain- wall vacuum solution [l|, [23 . The simplest version of the D=4, N= 1 action 
contains three chiral superfields, namely the dilaton 5, the universal T modulus and Z which specifies the position of 
the additional brane in the orbifold direction. In terms of the underlying component fields, these superfields can be 
written as |22| 

S* = e'^ + qz^e'^ - 2i{a - qz^x) , 

T = e/ + 2zx , (2) 

Z^e''z~2i{C-zx) , 

where if, /?, z, cr, x, d BuTC all real scalar fields, with (/?, /?, z having a geometrical interpretation in terms of the higher- 
dimensional theory. The field (p comes from the 5-dimensional dilaton, a part of the universal D=5 hypermultiplet, 
and measures the size of the internal Calabi-Yau three-fold averaged over the orbifold, so that the Calabi-Yau volume, 
with respect to a fixed reference volume v, is given by we"^. The dilatonic axion, a, will be the other zero mode 
surviving from the universal hypermultiplet. The /3 field corresponds to the zero mode of the (55)-component in the 
D ~ 5 metric and measures the size of the orbifold (more precisely the orbifold size will be irpe^ , where irp is a fixed 
reference size). Its associated axion, %, is related to the vector field in the 5-dimensional gravity multiplct. The field 
z originates from the world volume of the additional 3-brane and specifics the position of the brane in the orbifold 



direction; it is normalised as z G [0, 1], so that z = {z = 1) will correspond to the observable (hidden) orbifold fixed 
plane. The associated axion (, originates from the additional 3-brane and is related to the self-dual two-form of the 
underlying 5-brane. 

The Kahler potential which governs the dynamics of these three chiral superfields, has been computed in Refs |22.l2a | 
and is given by 

K=-ln(s + S-q ^^^^J ]-3\n(T + T) , (3) 



T + T 
at lowest order in e, the strong coupling expansion parameter 



e 



- n^P-V 



qe^-^ , (4) 



which measures the size of the string-loop corrections to the 4-dimensional effective action or the strength of the 
Kaluza-Klein excitations in the orbifold direction. From the 5-dimensional point of view we could understand it as 
the warping of the domain wall vacuum solution in the orbifold direction, or geometrically, the relative size of the 
orbifold and the Calabi-Yau. 

Considering the range of validity of our effective theory allows us to place constraints on regions of superfield (or, 
equivalently, component field) parameter space in which our solutions can be trusted. As stated above, the Kahler 
potential, eq.©, is computed to first order in the stong coupling parameter, e. Thus our model will only be valid in 
the weakly coupling regime, e ^ 1, where 

SH>^ + qTR {^^e'^>qe^). (5) 

In addition ensuring that a' corrections remain small, requires 

Tr>^ ( — /3>-ln(8^)) . (6) 

Finally as we have not modelled bulk brane - boundary brane collisions our effective theory will break down if the 
bulk 3-brane strikes a boundary, i.e. 

ZRe[0,TR] (^^ze[0,l]). (7) 

Subject to these considerations we now derive our effective theory. 
The four-dimensional N=l supergravity action, is of the form 

5 = ^ I x/^ Qi? + K,-d^<^^d'^^^ + V^ d^x , (8) 

where K.^j = ^.^^ is the Kahler metric (given in Appendix IX|) : <j>* = {S,T,Z) are the three complex chiral 
superfields; V($) is the scalar potential and Kp is the 4-dimensional Newton constant which, in terms of higher 
dimensional quantities, can be expressed as 

,.2 

SttGn , (9) 



2TTpV 

where k is the 11-dimensional Newton constant. 

Initial studies have been carried out upon the dynamics of this system in the absence of a potential. The full 
4-dimensional solutions were found in Ref. [5] and some aspects of the 5-dimensional solutions were discussed in 
Ref. [4]. In this article we will go beyond the perturbative approach, where V^ = 0, and consider non-perturbative 
effects, which induce an effective scalar potential whose most general form for 4-dimensional A^ = 1 supergravity is 



V = e^{K'W,WDjW-3WW) . (10) 

Here K'^^ is the inverse Kahler metric and DiW = diW + §§tW is the Kahler covariant derivative acting on the 
superpotential. 

There are a number of non-perturbative effects we could consider |25l \2dL 1271 123 . l29l l30| , each contributing to the 
superpotential for the chiral superfields. Possible effects include [2g: gaugino condensation, M2-branc instantons 
stretching between the two boundaries, an M5-brane wrapping the whole Calabi-Yau and M2 instantons that stretch 



between the nine-boundaries and the five-branes in the higher dimensional theory. This last contribution is amongst 
the dominant, therefore we will study its effect in our theory while neg lecting the others. Specifically, in our model 
with just one additional 3-brane, the superpotential takes the form J23 

WM2,A/5 = Me-'^^^ + e2^^(^-^^) , (11) 

where h is of the order of the Planck mass and, for practical purposes, we will use h — 1. The first term of eq. I|ll|) 
will correspond to the membrane stretching between the visible boundary and the bulk 3-brane, while the second one 
represents the M2 stretching from the bulk 3-brane to the hidden orbifold fixed plane. 

We can now compute the scalar potential, eq. (|10|l . for this particular case, which can be written as 

V^Vo{<PR)[Ai(t>R)+Bi(t)R)cos0{(t)i)] , (12) 

where 4ir={Sr, Tr, Zr) are the real parts of the chiral superfields, whereas (j}i—{Si^ T/, Zj) represent the imaginary 
parts. It is worth noting that only two of the imaginary parts of the superfields appear in the cosine term, 

9{Ti, Zi) = 2'Kq{2Zi - Ti) ^^ eix, C, z) = 2^q(-4C + 4zx - 2x) , (13) 

whereas the other terms in the potential depend solely on the real part of the superfields in a polynomial and 
exponential way. The full potential in terms of both superfields and component fields is given in Appendix IbI 

The action described in eq. © has the following symmetries 

Boundary- exchange symmetry 

The full action is invariant under following variations in the superfields, 

S — y S + q{T-2Z) , 

T — > T , (14) 

Z — > -Z + T . 

We can see from the previous equation that Z = T /2 is a fixed point under this symmetry, and this fact will be 
exploited when making our analytic approximations in section The geometrical interpretation of this symmetry 
becomes evident if we recast it in terms of the component fields. Inserting eq. ^ into eq. H14|l we find 

(15) 

and we see that our symmetry exchanges boundary branes. 

The Kahler potential, eq. Q, will always be invariant under this transformation however the superpotential, eq. Ijllll . 
is only invariant as we have chosen ft,, the prefactor in front of the exponentials, to be the same for both membrane 
instantons. 

Axionic symmetries 

When we write the action, eq. ^, in terms of the component fields we find, in the absence of a potential, two 
trivial shift symmetries in the axions 
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■^V ■ 


> CT - 


-> (T — / 


(i- 


->/? : 


, X- 


-^x , 


(1- 


--) : 


, c- 


^-c , 



cr, X~^Xi C-> C + ci , (16) 



and 



CT^CT + C2, X^X, C^C , (17) 

and a third, less obvious, symmetry which mixes the axions with z, the bulk brane position 

(18) 

This can be understood as the reduction to 4 dimensions ^ of the well known (e.g. Ref. |22j| ) gauge symmetries of the 
effective D — 11 action for Horava-Wittcn theory coupled to a five brane. 



a - 


-^ a + qc^z'^ 


X - 


-^ X + C3, 


c - 


^ C + C3Z . 



^ The associated 5 dimensional symmetries arc given in Ref. |22 



We saw in eqs. (|12() and H13() . that only two axions appear in the potential and that they are contained in its cos 9 
term. Therefore, although 8 is not invariant under the first, 9 —^ 9 + ci, or third, 9 — > 9 + 0^, symmetries, the potential's 
dependence on cos 9 means that we only break the full symmetries to discrete ones. In other words, ci = ±n/q in 
eq. p6|l whereas C3 = ±n/2q in eq. (|18|) . The second symmetry, eq. (|17|) . is unaffected as V is independent of a. 

III. NEW POTENTIAL FREE COSMOLOGICAL SOLUTIONS 

In Ref. |ll| it was shown that, when considering solely the kinetic terms, one can consistently truncate the component 
field action by setting a and C to constants and x to zero to leave the maximally truncated moving brane action. We 
can now use the third axion symmetry, eq. H18() . to generate a new, more general, class of analytic solutions, which 
include non trivial dynamics for the axions C and a, from the solutions found in Ref. [Tj. 

Before we generate our new solutions let us briefly discuss the solutions found in Ref. [il|. These will be relevant 
both in the context of our new axionic solutions and in certain limits of the potential driven system. They represent 
the solutions to the M-theory 4-dimcnsional effective action, explicitly given by 



^ ^ " 2^ / ^^^ (f + i^'^^^'^ + l^.Pd'^P + le^'^d.zd^^ 



The cosmological solutions in flat FRW space were found to be 

a = -ln\t\+ao , 



z]d^x . (19) 



p = p0,,ln|f| + (m/-mOln(|tr' + l) '+/3o , (20) 

_i 
LP = p^_iln|t| + (p^pj -p^,j)lnf |iP* + ij ' +ipo , 



z = d(l + |t|'') +Z0 . 
The initial and final expansion powers of the fields P and f are subject to the constraint 



3p|,„+P^,„ = 3> n^ij , (21) 

which defines an ellipse in the {pp,Pip^) plane, and are related by the linear map 

This ellipse and linear map can be seen in figure [21 where we have used them to understand the evolution of the fields 
in our potential driven solutions. 

The power i5 is given hy S = pp^i — p^p^i and can be fixed such that 

5 > (— ) timebranch , 
5 < (+) timebranch . 

Note that the initial position of the bulk brane, zq, is unconstrained, however the distance it moves, d, is constrained 
by the condition 

(po-/3o = ln(^^) . (23) 

In the asymptotic limits (t -^ ±0, t — > ±00), z is a constant and the fields ip and (3 behave as freely rolling radii, with 
their early (late) time expansion parameters given by ptp^i and pp^i (p^j and ppj) ^. In the intervening period the 



^ The parameters in both limits are constrained by eq. II21I however, at early times, they can only bo on the part of the ellipse I21i where 
Pf3,i ~Pv,i > Oi while at late times the fields satisfy ppj —Pt^.f < 0. Those negative time branch constraints get reversed for the positive 
time branch. 



3-brane moves significantly, leading to a more complicated evolution of the fields, but ultimately generates a transition 
from the early to late time rolling radii regimes. It is worth noting that the bulk brane does not necessarily collide 
with a boundary, it starts at zq and moves a fixed distance rf, to become constant again at zq + d. Thus, provided 
Zq + d G [0, 1], there will be no collision with a boundary. 

Now when we generate the axionic solutions from our new symmetry eq. H18|l . we see that their dependence on 
the moving brane ensures that the axions remain constant in the asymptotic limits, with the bulk brane movement 
generating a transition between these positions. Applying the non trivial symmetry, eq. (|18|) . to the solution for z, 
eq. I|2U[I. we see that the behaviour of the axions is explicitly given by 

a = qc^d^ [l + \tf^ + qc^zod (l + \tf^ +<Jo, 

X = C3 , (24) 

C = C3d(l + |i|')~' + Co. 

It is worth noting that, as in the case of the minimal solutions, these new ones have similarities to solutions found 
in the PBB literature. In Ref [Ij it was noted that the minimal moving brane solutions presented above, eqs. H20|l . 
are similar to the dilaton-moduli-axion solutions [33 of PBB cosmology |33, with the moving brane playing the 
role of the axion. Our new solutions, eqs. H24|l . in combination with the unaffected minimal solutions, eqs. 120|) . are 
themselves similar to the dilaton-moduli-'3 axion' solutions derived from an SL(3,R) invariant action in Ref. l37J . This 
suggests that, as in the PBB case, we would expect that our axions will provide us with a scale invariant isocurvature 
perturbation. 

In summary, although we generate non-trivial axion dynamics, the behaviour of the scale factor and the other, 
non-axionic, fields remains unchanged in this more general class of solutions. What we see is that the moving brane 
not only causes a transition in the moduli fields, as described above, but also causes the axions to evolve smoothly 
from one constant value to another. It should be noted that, as in the axion free case, the strong coupling parameter 
e grows asymptotically and the effective theory breaks down. At this point we can no longer say whether z or the 
axions remain constant. 

IV. INCLUDING NON-PERTURBATIVE EFFECTS 

We now return to the main topic of this paper, i.e. the inclusion of a non-perturbative superpotential into the 
dynamics of the 4-dimensional effective theory of heterotic M-theory. As discussed above we will only consider the 
effect of two membrane instantons, one stretching from each boundary to a single (moving) bulk brane. In this 
section we will concentrate on the component fields and their potential, as well as on simplifying the action as much 
as possible. 

From eq. (|12|l it is clear that all dependence on the axions is contained in cos{9{<f>i)) and, therefore, we will have 
a minimum of the potential along the imaginary directions for either 9 — or 9 ^ n. So, in principle, one could 
integrate out the imaginary fields and work with a simplified action. However, given that Vo{4>n)B{(j)ii) is a time 
dependent function of the real fields (p, (3 and z (or, equivalently, Sn, Tr and Zi{), we need to ensure that, if the 
axions are initially set to a minimum of the potential, this point remains a minimum. That is, we need to ensure that 
VqB, has the same sign throughout the evolution of our system. This is equivalent, see Appendix IbI to ensuring that 
the following function does not change sign 

M = {i- 2A'K'^qe^^f + IQ-k"^ q^ e^l^ z{z - 1) + Uirqe^) . (25) 

Independently of truncating the axions the validity of our effective theory already restricts the available (0, /3) space, 
as shown in eqs. ©-O. Then to truncate the axions our two possible choices are 

1 - M < 0. In the weak coupling regime, this condition is satisfied provided 

This leaves almost all of the remaining field space in which we can consistently truncate the axions. 

2 - A/ > 0. Reversing our choice leaves a very small region of (ys, j3 space in which we can consistently truncate the 
axions within a reliable 4-dimensional effective theory. 

Therefore, in the remainder of this section we will concentrate on the first, larger area of parameter space in which 
we can safely set x = 0; C = (or ±n/4g) and a = const. Explicitly our action now reads 

^ = " 2^ / ^^^^ (f + i^'^^^''^ + \d^.Pd^P + \eP-''d^zd>^z + F (^, /?, z)\ d^x , (27) 



where V = Voi<l)i)[A{(f)^) + B{4'i)], with;/), = ((^,/?,z). 

A. Equations of Motion 

We now look at the cosmological solutions of the system governed by our truncated action eq. H27|) . We consider 
the simple case of homogeneous, time-dependent, fields in a spatial flat FRW space time. Explicitly our ansatz reads 

0i = 4>iit) ■ 

As usual our solutions will have two branches, a negative branch with t, d < 0, and a future curvature singularity 
at t = 0_, and a positive branch with t, d > 0, and a past curvature singularity at i = 0+. In this section we will 
concentrate on the negative branch while in section when we include a perfect fluid, we will look at the positive 
one. 

In general the equations derived from our action cannot be solved analytically, and we will resort to numerical 
analysis. However we can sometimes make useful approximations. Inspecting the form of the potential given by 
eq. (jB4|l . we expect the dominant contributions to come from the double exponential in the (3 direction, given that 
the (fi and z directions are only proportional to single exponentials. Thus the potential acts to expand the orbifold. 
This suggest that, unless the kinetic energy or friction terms prevent the fields from running down the potential, the 
system will evolve towards the ^ = solutions (see section llTT|) at late times. 

In this context it is useful to parameterise the fields as RR fields, i.e. let 

if=p^\n \t\ + ipo , 
P = Pl3\n\t\+l3o , (28) 

with the RR parameters given by the expansion coefficients pp and p^. Using these equations we can retrieve effective 
RR parameters from the numerical solutions to describe the evolution of our fields. Recall that, in the potential free 
case, the fields make a smooth transition from one RR solution to another (defined by the map in eq. 1)22(1 '). which 
means that, although the expansion parameters are no longer constant when the bulk brane is moving, they evolve in 
a well-defined way from one value to another. Thus, provided the potential does not totally dominate the dynamics, 
the RR parameterisation of our numerical results will be a useful way to get a handle on the evolution of the moduli 
fields. 

It will also be useful to consider the shape of our potential in the regimes where the RR parameterisation is valid. 
When we substitute eq. H28|) into the potential eq. 1IJ4() . neglect the constants and note that z e [0, 1] in the regions 
where our theory is valid, we find that all the exponential terms have the same shape. Thus the form of our potential 
can be approximated by 

V{ip,P,z)K[{a\t\-^P0-P^ +b\t\-P0-P^ +c\t\-^P^')e-'^^*^''^] , (29) 

where a, 6, c, d G R^. 

For our numerical analysis it proves convenient to use a, rather than t, as our time parameter, although we will 
still use the time derivative of the fields to define our velocities. Practically this means we will numerically integrate 
the following equations 

'f' = -3^ - ^qe^-'^P - Iv,^ , 
a a 

/?' = -3^+^qe^-^z' - ^V,p , (30) 

6a Aa 

z' = -4(3d + /3-^)i-^e'^-'^F,^ , 
a qa 

where '= ^i ' — ^ ^-^d ^'U variables are functions of a. In addition we find that the Friedmann equation gives 



2- ^(1^2 + ^/j2^ 9/3-^.2^^)^ (3^) 



1.1 

3'4" 4' 2 



with a = — vq^ as we are considering the negative time branch. 

Within this framework we will not solve for the scale factor. This means that in order to retrieve the RR parameters 
from the numerics, we need to parameterise a = Paln \t\ + ag. We can then use the Friedmann equation to find p^ 
{Pa = d/d'). Note that p^ — > (1/3) when the potential becomes negligible and the scale factor has the usual (1/3) 
power law evolution of a system governed predominantly by kinetic energy, i.e eq. (|20|l . 



B. An explicit example 

To describe the effects of the potential on the dynamics of the system we present and discuss an exphcit example. 
In order to highlight most of the key effects of including a potential we have chosen an example where both potential 
and kinetic terms are important. We will then comment on some features not found in this example to develop a 
more complete picture of the model. 

In our example the potential plays an important role at early times without totally dominating the dynamics. Then, 
as it becomes less important, the system behaves similarly to the potential free case but punctuated by periods of 
more complex evolution, as the bulk brane approaches the boundary and is repelled by the potential. 

We start with a contracting Calabi-Yau (C-Y) , an approximately static orbifold and a slowly moving bulk brane at 
a = 3. Specifically we have 






3.1, (iin = 1.0, Zin = 0.3 , 

0.001, ft„ = 0.4 X 10"^, Zr 



0.6 X 10" 



(32) 



In the absence of a potential these initial conditions would lead to a very similar evolution to the explicit example 
given in Ref. [iJl, i.e. the C-Y continues to contract while the orbifold and bulk brane remain almost static until 
a « 0. At this point the bulk brane moves appreciably to its new stable value causing the other fields (}p and /3) to 
make the transition to their late time RR values. 

In this case the non negligible initial potential means that the fields do not behave as RR fields at early times. 
Instead the fields roll rapidly down the potential and go straight into the transition period where the bulk brane is 
moving appreciably and ip and /3 arc evolving towards their late time solutions. 

After this initial, brief period, the potential is only important when the bulk brane approaches the boundaries, and 
the fields evolve as in the potential free case. This is illustrated in figure ^ where we see that there is no early time 




FIG. 1: Evolution of the fields </?, /3 and z for the exphcit example given in the text. Recall that z G [0, 1] (dashed lines in 
plot) and that the negative time branch evolution means plots read right to left. 



RR evolution. The bulk brane is moving appreciably and the fields (^ and /3 start the transition to their late time 
RR solutions straightaway. From this point we see ^p and /3 continue their transition to the late time solutions while 
z rebounds across the orbifold, prevented from striking the boundary by the membrane potential, until it stabilizes. 

We can get a different perspective on this behaviour by considering the effective RR parameters. Figure|2shows the 
ellipse, eq. H21() . representing early and late time potential free solutions. As expected, in the presence of a potential 
we start well away from the ellipse, but as y — *■ 0, p^ and p^ return to its centre. They then make the transition to 
the late time half of the ellipse, with their track punctuated by four deviations as the potential becomes strong when 
the bulk brane approaches a boundary. 

Note that, although our numerical results gencrically show this rebounding behaviour, our effective theory docs not 
model a collision and is no longer reliable if the bulk brane strikes a boundary, rather than just approaching it and 
rebounding. 




FIG. 2: Evolution of the effective RR parameters for the example given in the text (solid line). The arrows map the negative 
time branch evolution between RR solutions in the potential free case, see eq. 1221 and the early (late) time part of the RR 
ellipse, eq. (I21II . are plotted with a dotted (dashed) curve. 



C. Interpretation of results 

In all cases we did indeed find that, at late times (i.e. i ^ 0_, a ^ — oo ), ip and (3 followed one of the y = late 
time solutions. Perhaps surprisingly, considering the shape of the potential, this was not necessarily accompanied by 
late time expanding orbifolds and the potential going to zero. 

In cases where the orbifold did contract at late times (parametcrised by p^ > 0) it was not free to contract as 
rapidly as in the potential free case. Using the RR approximation, eq. (|28|l . in the equations of motion ^ we find 
that the potential is negligible (even when growing asymptotically) at late times provided pp < (1/3). If the orbifold 
contracts any faster the potential acts to slow the contraction until pp < (1/3). This provides an extra constraint to 
our late time solutions on top of the V^ = late time conditions, i.e. and (3 will evolve to the part of the ellipse, 
eq. (|23, where 5 = pp — p^ <Q and pp < (1/3). 

These late time contracting orbifolds are of interest because the potential is no longer exponentially suppressed at 
late time and can grow asymptotically when pp > 0, see eq. H29|l . This could provide a way of transferring isocurvature 
perturbations into curvature perturbations in the moving brane scenario, a simplified version of this scenario is studied 
inRef. [s^. 

Finally, regarding our truncation, one should beware that, in the cases where p^ > 0, it is possible that, while in the 
weak coupling regime, we can evolve into the region of parameter space where the axions are no longer in a minimum 
of the potential, forcing us to include all the axionic equations. 

V. NON-PERTURBATIVE EFFECTS AND BACKGROUND FLUID 

In order to let the fields evolve in a cosmologically realistic background, we will include a background perfect fluid 
in the positive time branch of our model. 

We will now look at the solutions of the system governed by eq. (jH)) , with the particular scalar potential given in 
eq. ((T^ . It proves convenient to use the chiral superfield notation for the analysis, and then translate our results to 
the component field expressions. 

We consider the simple case of homogeneous, time-dependent, fields in a spatially flat FRW space time background, 
and let them evolve. Given this ansatz, we find the following equations of motion for the complex superfields 



$* + 3iJ$' 



r}j^$^$*= 



K'^djV 







(33) 



It is simpler to use the standard eom in terms of proper time derivatives rather than eqs. 1311 . 
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where $* = {S,T, Z), <j>* = d^^dt, djV = dV/dW and the connection on the Kiihler manifold has the form 

Tl=lC''—i^. (34) 

In addition, we obtain the Friedman equation for the Hubble factor H = a/ a = a, where a{t) ~ ln|a(i)| was defined 
in section ITVI and a{t) is the scale factor of the Universe, 

3i/2 = 4(p^ + Pb)= K,j^'^' +V + nlpb , (35) 

with pb the energy density of the fluid, which evolves with the scale factor in the usual way 

Pb = P6o/a'<"^+'^ , (36) 

where Wb defines its equation of state: pb — WbPb- We will use radiation, i.e. Wb = 1/3. 

It is worth splitting the equations of motion for the complex chiral superfields into real and imaginary parts 

4 + 3ff0^ + \t)S'r^r - 0i07) + \K'^d,,V = , (37) 

<^} + 3H^} + ir;., {iP.^j,), + ^=^4>1) + \K'^d,, F - , (38) 

where now (/)^ = {Sr^Tr^Zr) {(j)\ = {Si,Tj,Zj)) refers to the real (imaginary) part of the superfields; djj^ [dji) are 
used to denote the derivative of the potential with respect to the real (imaginary) parts of the fields. The scalar 
potential V is given, in terms of chiral superfields, in eq. (|B2p . As already mentioned at the beginning of section IIVI 
the dependence of this potential on the axions is contained in the cosine term, and a minimum is found for 9 = ±2n7r, 
i.e. 

Zi = \{Ti±^) , (39) 

z q 

and no explicit dependence on 5"/. It seems natural, therefore, to start by solving eq. (|37() with the axions fixed at 
the minimum shown above. 

A. Truncated action 

We start off by choosing a particular vacuum for the axions, which will be given by 

(40) 



Si-- 


- Sio 


Ti-- 


= , 


Zi 


2q 



Note that the real parts of the fields are bounded by the constraints that define the field-space region in which the 
low energy effective theory is valid, written in eqs. ®-([7I). We restrict our model to the M < region, where M was 
defined in eq. lP5|l . to consistently truncate the action. 

We proceed by evolving the real parts of the superfields in the presence of a background fluid that dominates the 
total energy from the beginning. Under these initial conditions, the fields undergo two different stages, as we can see 
in figure |31 

1. At early times the fields start evolving because of the initial velocity they are given, which is much more 
important than the potential energy they have, until the friction term 3H(j)' becomes big enough to freeze them. That 
is, friction holds the real parts of the superfields on the potential slope where they remain fixed at a constant value 
for a while. 

2. As the background fluid energy density, pb, decreases with time, the Hubble factor becomes smaller, and the 
friction term is no longer able to hold the fields. The real parts of the superfields then roll down the potential slope 
and are driven to an attractor which makes them scale like the background fluid, as we can see in figure 21 Although 
the fields do have an attractor solution which makes them mimic the fluid behaviour independently of their initial 
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FIG. 3: The real parts of the superfields, Sr, Tr, Zr as a function of a{t), which is essentially the number of e-folds. The 
dash-dotted lines show the analytic approximations explained in the text. 




5 10 15 20 25 30 35 



FIG. 4: Logarithm of the energy densities for the background fluid, pt, and the real parts of the superfields, p^, as a function 
of a. 



conditions, they never dominate the background energy density. This, as shown in Ref. 34\, rules out any hope for 
them to become the quintessence field. 

In figure |5l we plot the equation of state of the real parts of the superfields as a function of the number of e-folds, 
showing that, at late times, they will oscillate around w ^ Wh {= 1/3 throughout this paper). This, again, is far from 
the quintessential behaviour, where —1 < wq < —0.7. 

We can even go further and provide analytic solutions for these two stages. 

1. Early stage: the evolution of the fields is dominated by the kinetic term, thus we can neglect the scalar potential 
(V ^ 0). In this regime we can use the complete solution for the fields and their equations of motion given in Ref. |38l | , 
but under the assumption that the background fluid dominates in the Friedman equation, i.e. Pkin / P fluid <C 1. The 
early behaviour of the fields can then be approximated by 



n 



Ai +A2e- 
Bi + Bae" 



(41) 



Z^ « mz, + Z2e-n • 
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FIG. 5: Equation of state for the real parts of the superfields as a function of a. 



Therefore, as a grows the second term on the right hand side of these equations decreases quite fast, and we end up 
with the constant values that correspond to the freezing period we mentioned before. 
2. Late stage: in the scahng regime the evolution of the fields can be approximated by 



(42) 



It is also possible to determine the value of these six constants. To do so, we will perform some approximations in 
the equations of motion for the real parts of the fields, eq. (|37|l . which, in terms of q, is given by 



s^^. 


. Cia + S^^ , 


Tk^^ 


. C^a + T^ , 


zk^ 


« C^a + Z™ 



where </>' = d4>/da. These are the following 

• The background fluid dominates the total energy density, which results in two conditions 
a) Its energy density, pf,, is much bigger than the kinetic energy of the real parts of the superfields 



(43) 



b) It dominates over the scalar potential 



V < n^ppb 



Therefore, the fluid will dominate the Friedman equation 

3i?^ 



2 
KpPb 



(44) 



(45) 



(46) 



The symmetry in eq. (I14II for the real parts of the superfields, has a fixed point at Zji = TbJ2 which defines 
a minimum of the potential along the Zp direction. In the positive time branch case, the presence of the fluid 
acts as a friction term that causes the fields to have a kinetic energy of the same order of magnitude as the 
potential one. This is why the Zr modulus always ends up along the valley that defines the minimum in the Zr 
direction. Physically, this means that the additional bulk 3-brane will always be stabilized at the middle point 
in the orbifold direction, i.e. in the middle between the two boundaries (observable and hidden). In this sense, 
the background fluid stops the fields from acquiring enough kinetic energy to jump out of the valley defined by 
minimum of the potential in that direction. 

We can therefore reduce our system to only two real fields by substituting Z|j ~ 2-^i? ^^^ '^^^ equations. 



13 

• We arc in weak coupling regime (e ^ 1), so that ^Tr ^ Sr. 

• The mixing terms are small for the region we are working in, i.e. T" 0}j0]^ <^ 30^ . 

• The slope of the potential in the T/j-direction is much steeper than in the S'j^-direction as the potential depends 
exponentially on Tj^ while it has a polynomial dependence on Sn, therefore K^'^"dTf,V ^ K^^djV for j ^ Tr. 

Equation (|43|l will, after making these approximations, reduce to 



m^{i-wb) 



(47) 



=Tr/2 



From this we can sec that the different scaling slopes for the different real parts of the superfields will just be due to 
the K^^'^ term, which will determine the following relations 



K 



SrTr 



K 



ZrTr 



= 


_ 1 j^TnTn 
4 


= 


=> Ci = 


-Ic, 


ZR=Tn/2 


Zn=TR/2 




Zr=Tr/2 


2 


Zr=Tr/2 


=> C3-- 


-¥^ 



Now, solving eq. H47(l together with eq. H46|l . we can easily obtain 



Cl = 


^(l + .«.) 


Cz = 


2',^^ + -) 


Cz = 


k-^^ 



(48) 
(49) 

(50) 



and 



Tiri 
c 



1 



2^2 „3 



In 



Mh^Ti^q 



27Tq 27(1 - wl)poKJ 



(51) 



Note that this expression depends on q, tob and pbo, so that Tr will end up having the same values independently of 
the initial conditions for the moduli. This implies that the orbifold size will be, at any time, dependent only on the 
characteristics of the background fluid and on the charge of the moving brane. 

We can now match both stages: using the approximations in both regimes and matching them for a 3> 1, we can 
compute the rest of the constants in eqs. (|41|) and (|42|) . and have a complete analytic solution given in Appendix lUl 
together with an estimate of the time at which the scaling regime starts 



a 



1 



(Bi ~ Tl 



To complete the late stage approximation, we will write also the matching constant for the S'jj-field 






S^,^^{T^-B,) 



(52) 



(53) 



where all the remaining constants are given in Appendix lO 

As we can see in figure|21 the analytic solution reproduces accurately the numerical one for the early and late times. 
There is a slight discrepancy at late time in our analytic approximation for the S^, but we can see that the slope of 
its late time approximation, Ci, is correct. 

We can now express these solutions in terms of the component fields; for the early stage of the evolution we get 



Ze. 






(54) 



whereas, at late stages 



^Pl - 


» Ci , 


A ^ 


» Ina 


Zl ? 


1 
"" 2 ' 
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C2 



(55) 



where all the constants are given in Appendix [CI In figure we have plotted the evolution of the component fields 
together with the early and late time approximations. Again the analytical results coincide with the numerical ones. 
We can, for completeness, study different initial conditions to the background domination scenario assumed until now. 




10 15 20 25 30 35 



FIG. 6: Evolution and analytic approximations (dotted lines) for the component fields (p, /3, ^ as a function of a. 



1. Starting with a moduli dominating energy density 

This will translate into the condition, 

PbO < PipO 



(56) 



constrained to px = Pb + P4> < Mq, where Mg defines the GUT scale. The real parts of the superfields will, in this 
case, evolve quite fast, but the friction term due to the fluid will slow them down so that, at some point, the field 
energy density will get smaller than the background one. Thus the fields will either grow too fast to stay in the weak 
coupling regime, or follow the behaviour explained before, i.e. once their density becomes smaller than pt they will 
end up scaling with the background energy density. The only difference will now be that we cannot use the early 
approximations anymore. 



2. Starting with a large kinetic energy for Zj^ 

When we give the Zn-iield a large enough initial velocity, the 3-brane will move towards the boundary and hit it, 
bouncing back again, causing our theory to break down. There are certain cases in which the 3-brane will just get 
close to the boundary, go up the potential slope and bounce back. This will affect the evolution of the other real 
parts of the superfields because of their mixing terms in the equations of motion. In the case of Sr, which is coupled 
to Zn in quite a strong way through T^" z ^^ ^"^^^ ^^ ^s"z ' ^^^*^ ^^^^ ^^^^ ^R bounces back will drag it along, so 
that the volume of the Calabi-Yau will become smaller. Meanwhile the Tfj-modulus will not be affected by the other 
fields' behaviour, as in its equation of motion the dominant mixing terms is by far T^^^, . The decrease in the value 
of Sr, while Tfj is still increasing, will lead the system to strong coupling (that is Sr 3> qTji implying e S> 1), so that 
our model breaks down and results are not valid anymore. 
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B. Including axions 

If we let the axions evolve as well, we will see that their behaviour is similar to that of the real parts: they will 
just evolve while their kinetic energy is bigger than the friction term due to the background fluid and then, they will 
slow down and be driven towards their minimum where they will stabilize at their minima after some oscillations, as 
we can see in figure The evolution of the real parts of the superfields will not be affected by the axionic one as the 
axions only appear in the real parts of the superfields equations of motion, eq. (|37|l . through the mixing term, which 
is typically very small. Only in certain cases, in which the axions have a large amount of initial kinetic energy, will 
the real parts of the superfields modify their behaviour. Because this mixing term has a negative sign, the real parts 
will decrease while the axions still have an effect on the evolution, until the friction term dominates again pushing 
the axions to their minimum, and letting the fields go to the same scaling regime again. As before, the slopes of this 
late approximation are really accurate. 




FIG. 7: Evolution of the axions, Si, Ti, Zi as a function of the number of e- folds, a. 



C. Interpretation of results 



There are four main results we would like to highlight: 

1. The axions, in general, will not affect the real parts of the superfields evolution, as they will always stabilize at 
the minimum of the potential in the axionic direction, thanks to the damping of the background fluid. Thus, when we 
consistently truncate the action to its minimal version, where the axions rest at their minimum, with certain constant 
values given in eq. H4(J|) . we will not lose any important features. Only if their initial kinetic energy is large will they 
have an effect on the real parts of the superflelds evolution, but just through the early stages, leading to the same 
scaling solution at late time. 

2. The additional bulk 3-brane gets stabilized at the middle point between the two boundaries, which is shown in 
figure El by the z field going always to the value z = 0.5. This behaviour is due to the symmetries in our potential, 
eq. (|14() , together with the fact that there is a background fiuid that dominates the total energy, preventing the kinetic 
energy of the fields from growing over the potential one. 

3. The orbifold size, TrpTu, grows linearly with the number of e-folds, eq. H42|l . with a speed, C2, that depends only 
on the charge of the moving brane and the background fluid properties, as shown in eq. H50() . Thus, its size will be 
independent of the initial conditions of the chiral superfields. 

4. Our model always tends to strong coupling regime at late times. This was suggested in earlier papers in which 
the moving 3-brane was free, as there was neither a background fluid nor an instanton superpotential, see Rcf. [ll|, 
so we can conclude that the addition of the membrane instantons does not affect that result. We can, in this case, 
understand why this happens: it is due to the slope of the T^-modulus (C2) being bigger than the one of the S'^-field 
(Ci), which comes from the fact that X^«^« is bigger than K^"-'^" in this regime where Zfi = 7^/2. 
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5. The volume of the Calabi-Yau. (j), tends to a constant at late time, as we can see in eq. H55|) . This is due to the 
relation between the scaling regime slopes Ci = (7r/4)C2, that will cancel the time (a) dependence in (/>; 

</); = In (5j, - gZJ^VrA) = In {S'^^ |r™) , (57) 

but this will happen when going to strong coupling, where our theory breaks down. 

VI. CONCLUSIONS 

In this paper we have explored further aspects of low energy heterotic M-theory. We have presented new symmetries 
of the full action and used these to generate a new class of solutions with non-trivial axion dynamics. We saw that the 
non-axionic fields have the same evolution while two of the axions now make a similar transition to the bulk brane, 
from one constant value to another, leaving the system to return to a two field rolling radii model asymptotically. 
The similarity of these solutions to certain PBB ones suggests the axions will produce a scale invariant isocurvature 
perturbation. 

We then investigated the effect of membrane instantons on the dynamics of the system, both with and without 
a background fluid. In the absence of a background fluid we find that, despite a more complicated evolution, the 
Calabi-Yau and orbifold moduli behave as rolling radii in the late time limit with the potential dropping out of the 
dynamics. During the evolution to these late time solutions we find the membrane potential attempts to prevent 
the bulk brane striking a boundary, leaving it to stabilize in the bulk. Additionally, although negligible in the 
background equations of motion, the potential can grow asymptotically at late times, providing a possible mechanism 
for transferring isocurvature perturbations into curvatiu:e perturbations. 

We then included a perfect fluid in the background for the positive time branch of the evolution of the Universe. 
This fluid acts as a friction term for the fields evolution, stopping their kinetic energy from growing bigger than the 
potential one. This has an important effect in 

i) The axions which, after evolving, stabilize at the minimum of the potential in their direction. They drop out of the 
real part of the fields equation of motion H37|) and, while evolving, do not change the results found for the truncated 
version of the action. 

ii) The moving brane, that is driven to the minimum of the potential in its direction, that is, it is stabilized at the 
middle point between the orbifold boundaries. 

In this case, the real parts of the supcrficlds scale with the background fiuid at late times. This is a nice feature, 
typical of any quintessential model but, in this case, it is insufficient to explain the amount of dark energy in the 
Universe, as the moduli energy density never dominates over the background fluid one. We were able to find a set 
of analytic solutions, given by eq. (|42|l for all the fields. In general Sfj, Tr and Z/j increase linearly with a, with a 
speed, eq. (|50|l . that depends only on the charge of the 3-brane, q, and the background fluid characteristics, pB, and 
which is independent of their initial conditions. 

Because of this scaling behaviour, the orbifold size, npTu, ends up being always the same at a certain time, given 
a flxed background fluid and charge, as we can see in eq. H51|l . But the different slopes of the scaling solutions for the 
different fields drive the system to strong coupling at late times. 

In the future, we will investigate two possible directions. Firstly, we will try to find the higher dimensional fields 
that will act as a background fluid in our 4-dimensional effective theory and see the changes this might have in the 
superfields evolution. Secondly, we will be calculating the perturbation spectra of two of the new classes of solutions, 
i.e. the new axionic solutions and the solutions with a potential but without a fluid, to see whether either can generate 
scale invariant perturbations. We know that the minimal, potential free, solutions [ij cannot generate either adiabatic 
or isocurvature perturbations with a scale invariant spectrum Ref. p32 |. 
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APPENDIX A: KAHLER MATRIX 



The Kahler matrix, K^j — „!^. £v^. , takes the specific form shown below when using the tree level Kahler potential, 



eq. ©, 



/ 1 



^% 



QZr 



qZ% 1 (q^Z% ,„oZ% 3,,2^ 



_2qZR 
Tr 



v-^ -#(2/+^) ^(?/+^) ; 



where y ~ 2SR — 2qjfi^. Note that this matrix is a function of the real parts of the fields only and that it is symmetric. 

APPENDIX B: SCALAR POTENTIAL 

The D=4, N=l, supergravity scalar potential, eq. (|10|l . for this particular Kahler matrix but for a general super- 
potential, is of the form 



V 



1 



8y3r3 



\yWs~W\'' + ^ 



ZrWs + ^Wz 
2q 



TflWr + ZrWz + ^Ws ~ -W 



2,\W\ 



(Bl) 



where y has been defined above and Wi = dW/d^'^. 

When we include the specific superpotential of eq. 1)11(1 . as it was already presented in the text, the scalar potential 
in terms of the chiral superfields can be written as 



where 



V = Vq{Sr, Tr, Zr)[A{Sr, Tr, Zr) + B{Sr, Tr, Zr) cosd{Ti, Zi)] 



(B2) 



Vo{Sr,Tr,Zr) - -_^^_^__^-^ , 

A{Sr,Tr,Zr) = (-3 - 24^ VflTi? + S^'g'^fl - 247rgZj,)e-4^«^« 

+ (-3 - 24TT^qSRTR + SiT^q^Zl - IGtt^^TI 

+ 32Tr^q^TRZR-2ATTqTR + 2ATTqZR)e^'"i^^''~^''^ , 

B{Sr,Tr,Zr) = 2(-3 + 24^VflTfl - 8^V^^ - leTrV^flZ^ - 12^r7rfl)e-2-9^« 

eiTi,Zi) - 2Trq{2Zi-Ti) . 

Inserting the definition for the component fields given by eq. |(2J) we have 






48e3/3+0 ' 

(-3 - 2An'qe/+'^ - WTT^q^e^^iz - if + 24nqe''{z - 1)^-1^' (^-^'> 
+ (-3 - 24n\eP+^ - Wn^q^e^^iz)^ - 24Trqe^'z)e-^^'J^'' , 
B{(j),P,z) = {-6 + A87r'^qe/+'f' + 327T^q^e^^z{z~l)~247rqe'^)e-^'"^''^ , 
eixX.z) = 27rg(-4C + 4zx - 2x) • 



(B3) 



(B4) 
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APPENDIX C: ANALYTIC SOLUTIONS 

The different constants appearing in the early time approximations for Sr, Tr and Zji, eqs. H41|) . are given by 

Ai = Sr{0)^A2 , 

, snioy 

c 
Bi = Tn{0) - B2 , 

i>2 = , 

c 

^1 = ^ , (CI) 

BI + B2 

C = -{l-Wb) . 

The anafogous expressions in terms of the component fields 0, /3, 7, see eqs. H54|l are 

a, = ln{S^-q^)-a2 , 

(29 — ' ■ '■ r ■ — , 

^"^ [^r-Ir - q^ji ) 
61 = lnTir-62 , 
_ -Th(0)- 

°2 - ^rp,n ' (,C2j 

Z^ 

z^TBioy - Znioyr^ 

Z2 — T ■ 



cT"" 



Finally, the constants for the late time approximations in terms of the component fields, see eqs. ()55|) are 



4 
C2 - lnC2 . (C3) 
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